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ABSTRACT 


Wave spectra calculated using the Parameterized 
Nonlinear Wave Solution developed by Le Mehaute et al. 
(1984) are compared with field data acyuired at Leadbetter 
beach, Santa Barbara, California. The parameterized solu- 
tion satisfies the nonlinear free surface boundary ccndi- 
tions to a specified degree of accuracy and is expressed in 
terms of a converging truncated Fourier series. The wave- 
number, surface prcfile and wave orbital velocities are 
determined by the wave height and wave period at the local 
depth ci water. Spectral components are compared between 
the model results and field data. Good agreement is cbserved 
for waves corresponding to Ursell numbers ranging from 25 to 
75. #j;$For large Ursell numbers (strong nonlinear effects) the 


parameterized model underestimates the data. 
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The random nature of ocean waves is most often descrited 
uSing the linear spectral model, and many engineering frob- 
lems related to offshore and coastal structures are solved 
by employing the concept of directional wave spectra. Waves 
in the sea, however, can exhibit nonlinear properties. A 
departure of surface elevation from the Gaussian distriku- 
tion with positive skewness iS an example. Phase velocities 
of high freguency component waves not satisfying the linear 
Gispersicn relation is another example. Breaking of waves 
is a spectacular exazple of nonlinear benaviour of water 
waves. 

The study of nonlinear, water waves traveling over a 
horizcntal bottom dates back to Stokes(1847), who solved 
this froklem in the form of a power series in terms of a 
Stall parameter '€' related to the average wave slore (€ = 
ak, where ais the wave amplitude and k is the wavenumcer). 
it was fcund, however, that the Stokesian type series are 
honuniformly convergent, and are only valid in deep and 
intermediate water defth. 

Boussinesg (1877) and kKorteweg and de Vries (1895) 
develcped the cnoidal wave theories, which are based urfon 
power series in terms of wave height relative to water 
depth. These power series are uniformly convergent in 
Shallow water, but are invalid in deen water. 

Making the connection between the two above theories, 
Goda (1983) developed an empirical parameter to describe the 
phenotena of wave ncnlinearities with good resuits. The 
parameter bridges the wave steepness in deep water and the 
Ursell's parameter Ur = (a/d)/(d/L)2 (with a the wave 
amplitude, d the water depth, and L the local wave length) 
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in very shallow water, Spanning the full range of water 
waves. 

When ccnfronted with experimental results, Le Mehaute et 
al. (168) demonstrated that the nonlinear theories were not 
much rFetter than the linear wave theory; nonetheless there 
1s a need to clarify the extent of wave nonlinearity in the 
sea so that engineering problems can be solved with much 
more confidence and accuracy. 

During the last two decades, numerical approaches kased 
on truncated Fourier expansions have been proposed and well 
verified experimentally. The first such approach was devel- 
oped by Chappelear (1961) involving the use of the velccity 
potential. Dean (1965) used the stream function to develop 
his numerical wave tkeory, which was computationally simpler 
than Chafrpelear's technigue. Dean's stream functicn fcrn 
satisfies the Laplace equation, the kinematic free surface 
boundary condition, and the bottom boundary conditicn; the 
parameters in the stream function expression are chosen Ly a 
best fit to the dynagzic free surface boundary conditior. Von 
Schwinc and Kkeid (1972) developed characteristic solutions 
to finite amplitude waves using the velocity potential. 
Cokelet (1977) extended the method originally developed by 
Schwartz (1974) to allow a very accurate calculation of the 
characteristics of water waves. The procedure invclves 
expressing the complex potential solution in a Fourier 
series and represents the Fourier coefficients in terms of a 
perturbation parameter. Rienecker and Fenton (1980) used a 
finite Fourier series to describe waves by solving a Set joe 
nonliné€ar eyguations using Newton's method. Recently, Le 
Mehaute et al. (1984) developed a parameterized solution 
which is valid For all practical ranges of values of wave 
amplitude, frequency and water depth. In view of its rela- 
tive Simplicity, and experimental validity, the parameter- 


ized sclution iS reccmmended for engineering applications. 


The objective of this work is to utilize the spectral 
model Ey Le Mehaute et ale (1984) to compare with wave data 
acquired at Santa Barbara during the 2-6 February, 1980. 
The range of application of the model will be checked 
against measured wave spectral components and individual 
waves. The wave nonlinearity and Ursell number are analyzed 
uSing the model. Ccmments and conclusions are made for 


further study. 





A. PARAMETERIZED SOLUTION TO THE NONLINEAR WAVE PROBLEM 
1. Cverview 


The higher-order Stokian wave theories become alge- 
braically very complex. For practical purpose, it was 
desirable to develop wave theories that could be computed to 
any order. Le Mehaute et als (1984) developed a parameter- 
ized solution to the monochromatic, irrotational, nonlinear 
water wave over a acrizontal bottom. The formulation is 
expressed ina closed form in terms of a truncated Fourier 
series. All coefficients are expressed in terms of simple 
analytical functions which contain three parameters, namely, 
wave height, H, wave feriod, 1, and water depth, d. The 
number of terms in the Fourier series 1S proposed paramet- 
rically fora given level of accuracy. The solution is 
assumed to be applicable for all practical ranges from deep 
to Snallow water waves, i.e., d/Lo 2 0.005. The deep water 
Soluticn for d/Lo = 0.5 1s assumed to be valid for  d7iame 
0.5. 


As in the periodic water wave problem the formula- 
tion has iinear and nonlinear boundary conditions and a 
iinear geverning difierential eyuation. 

Chcosing a coordinate system moving with wave speed 
c= O/ k sotnat the time dependency is removed, i.e., a 
staticnary problem , it 1S possible to asSign a COn=aaue 
value to the stream function at the free surface. The vari- 
able Q is the wave angular frquency( 0 = 2 71 /T) and k is 
the wavenumber (k = Z7{/ L). The solution is a velccity 


potential function of the forn 
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D {(x-ct),z} Cr D (GO 2) which satisfies the Laplace 


equation : 


V7 7@,> 9 (2.1) 


The variable g =Soreegmeteaenotes phase, and the z axis is 
vertical positive upwards from the S.W.L. (still water 
level) and the horizcntal x axiS 1S positive in the direc- 
tion cf wave propagation. 


The kinematic bottom Ecundary condition 1s written 


as 
| = 0 DeayD 
Q, a (2.2) 
The dynamic free surface boundary conditicn, is 
derived using Bernoulli's equation written in the dimen- 


Sionless form at Z = i), 


2 


2 
= n+s-\($,-9) +b;|- 55 = ey = 10) 


VS (22) 
where H is the wave height, 7) is the water wave surface 
elevation and Be = an arbitrary Bernoulli constant. 


The kinematic free surface boundary condition is 


up # ceo = eo (OQ) (22) 


The e's are the allowable errors which snould be small. 

The measure cf how well these two boundary ccndi- 
tions are satisfied is defined by Ey. and E 5 , which are 
the mean sguared errcrs of the dynamic and kinematic free 


Surface Ltoundary conditions : 


1 


af 1/2 


oe ae é 
eS (8) (2.5) 


J ; 1/2 
poe aa: e 
z= —) 2; (0) (2.6) 
Jima 
in which j = sample at evenly spaced phase angles (€.g., 
with 5° interval) along the wave ;rofile. The Bernoulli 


constant, Be, 1s determined such that the variation cf the 
phase speed c is minizized over a wavelength. For an exact 
SOlUucion: E, and E, weuld be zero. 

A solution in the form of a limited series (or trun- 


cated Fourier expansion) is assumed 


O= - ae O y cosh [ak (d +2) | sin nQ 


(2.7) 


in which oP = dimensionless potential function; An = a 
derived coefficient; nis an integer and N = the number of 
terms. All the internal wave <field terms (velocity ccmzo- 
hents, dynamic pressure and acceleration components) can be 
Girectliy determined from this equation. EQ. Z./7 Satisives 
Eoth the Laplace equation (2. 1) and the kinematic botton 
boundary condrtroene( 2. rec aetui ys 

The equation fror the free surface correpending to 


60s. 2.7 JS 





N 
mia K ans >a sinh nk (d + 7))) €0s nQ 


(2.8) 


in which K = an integration constant satisfying 


T 
| mi = © 
: (2.38 
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At the wave crest ({ a See) = Tle and the wave trough 
(O=T7 ), 7) = 1), +Therefore 


We ‘ | ie i +E, 


7 (2% 10) 


on wick E 3 is to be determined and must Le small. 

Le Mehaute et al (1984) defines the phase velocity c 
or wavenumber k in accordance with the Hamiltonian varia- 
tional frinciple of rinimum energy, which, in practice, is 
oktained through minimizing £, with respect to the waven- 


umber, k 


dE, 


Ok = © (eal) 





Such a system of ecuations can eventually be solved asa 
function of phase velccity. 

The solution to the preblem consists in providing an 
analytical solution ‘which describes all the wave character- 
istics as functions cf oniy three parameters : wave height, 
H, wave period, T, and water depth,d. These parameters are 


grouped into two dimensionless parameters 








Z 
Dezth parameter; O= = en (21) 
Le rosea 6a | 
a 
H a0 
Wave steepness parameter: = = Mets 
T= em TT g — 
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in which Lo is the linear deep water wavelength and a = H/2 
is the wave amplitude. 

Having specified the initial boundary value prrecblen 
and analytical forms cf eqs. 2.3 and 2.4, the proklem now 
consists of expressing all unknowns ; An, k, K, and WN (eq. 
2.8), in terms of the two parameters 6 and 2 , such that 
a desired level of accuracy is achieved. These unknowns are 
determined by satisfying the boundary criteria as closely as 
possitle. Eq. 2.11, which specifies the phase velocity and 
wave number, must alse be satisfied. The method of solution 
is one cf trial and error involving multiple iterative 
Frocesses aS expressed in the flow chart in Figure 1. 

For practical pirposes, the final formulation nay 


be presented in the fclliowing: 


INPUT CATA : Given EF, T, and d then 


lox gT/2TT , MN arte kOe : C=2i7 a 


ae 


O=d/Le and y =H/L, _ Also. B= y [tannca. 5880 


a ~ 
ag O = gk tanh bigs. 


The symrkol * denotes the linear solution). 


FORM CFE SOLUTEOwsS: 


N 
waa oa [7k (d+2z)] sin nf (2. 14) 





2) 


N 
Sa usla(eenjorns 
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INTEGRATION CONSTANT: For application purposes, a parame~ 


terized form of K is determined by 


een 7G) ak onan 


Gee (2. 16) 


The use of K as expressed aktove generally yields tJlarger 
values cf E's than the value of K obtained by eg. TES, 
eicectly. 
NONLINEAR WAVE NUMBER : k = Z7T/L 

ee 0 Sen 0- 1 


b*® 


0 a Ep 
| | + 1879/3" °°" _ 9.09437 in i2 «10°51 Cae 


2) ee 
.e 


) 4 3466 qapee (2. 18) 
tanh k¥d 


k = po 





mene socluticn valid fcr 8 = 0.5 is also considered valid for 
b>0.5 provided that all the parameterized functions which 
follcw are determined with } = 0.5 irrespective of their 
real value; this means that deep water solutions ccrre- 


Sponding to 6 > 0.5 are not influenced by water depth. 
Seer FICIENISs 


Coefficient AY: 


A, = A™(1-R) 
(2. 19) 


“cosh Kg (2a 


} b, 
R= a (5.882 /? ) / 


(2.290 
ip which at = 0568 132" Sep 0c tSe [anc19° 6 ns 
bY = 0.0747 exp {1.1957 fin (10” 6 yo; 
2. © Ol1) < COmanONS 
4 
R = 0.5+(-03+¢))exp) 4 [101061] | oe 
PW wWhaGh a = ef aaa 
pu = 1.369 tanh (24.05/37 °° y 
c4’ = 1.3676 exp (179.037 °° 
Coetficient Az : 
ie O05) < iG Oe 
Ay = 0.4974, exp | a, (5.882 (3 )°2) 2 
ip when a‘, = -0.00633 exp {1.494 fan (10> ) yo"? 
bs, = -0.953 exp{-0.029 [ An (10° 6 "> 
oe COMI oO aD 
z 
p= 2.5x107 A, exp [a exp (b, 6°7)] (2.24) 


-0.98 
in which a Ias ey exp (~0.0134 /9 ) 
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: iS 
iE] — -_ —_ 
b", 16+15.318 exp ( 0.0017/3 ) 
-/.04 
= -~-U,. = e 
c%, 5-4.328 exp (-0 0112/7 ) 


Coefficient Az : 


feo-005 < 0 sano. t 


A, = 0.6724, exp a, (5.88273) °3) 


(2.25) 
°78 
in which at, = -0.021 exp{1.355 [2 (10° 6 i}” } 
wntze 
b', = -0.859 exp{-0.0747 [4n (to 6 ) } 
2 aio aS 
A = A exp (ay 
3 < 3) (2.26) 
0.52 .-0.4/8 
in which a" = Ba S50 B 
Boetrrcient An {nm > 3) 
/ j 7 US 
An = A = exp - 3. 5) 
ay ib ) (2a07) 


(ad 


- Analytical Validity 


The validity cf this theory is defined through the E 
Memmies as derived by eqs. 2.3, 2.4, 2.5, 2.6, 2.1 and 2.11 
Ey order of importance. The magnitude of the E's depends on 
the number of terms N used in the series. Beyond a certain 
value of N, it is found that the E values tend towards 
asymptotic values, indicating that they are not reduced by 


adding more terms, i.e€., the series is bounded. With this 


21 


Froperty in the solutions, the number of terms N reguired 


can be expressed in farameterized form by : 


N = integer part of{u + n, (5.38 B - 9.2)} (2423) 
n= 1+ 19.13 exp {-9.344 [1a (193 fe) oe (2azsh 
n,= 32.72 exp (-0-212 [in (105 () jit * a9 (2. 30) 


Le Mehaute et al, (1964) showed that the parameter- 
ized soluticn approxinates the free surface dynamic boundary 
condition as well as, if not better, than any other existirg 
theory. It also approximates the kinematic condition gener- 
ally Letter, except for the stream function theory devel 
oped Tv Dean (1974) by virtue of its definition. 

iruncated Fourier series do not represert lizit 
waves well since tne cusp at the wave crest reguires an 
infinite number of terms as was demonstrated by Dean (1974). 
Therefore, the accuracy of the parametric solution decreases 
rapidly as the wave heigkts approach limit wave conditions, 


which is where Ur is large. 


Be MEASURE OF NONLINEARITY 


The parameter which represents the wave nonlinearity in 
deep water is the wave steerness, or a measure of the 
Surface elevation sicre. The most common representaticn 1s 
either the ratio of wave height to wavelength H/L, or the 
Parameter ak. 

In the region of shallow water waves, the Stokian wave 
prozile for higher crders 1s an expansion using the nondi- 


mensicnal term 


Ur = (a/d)/(d2/L2) (2.31) 
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as the perturbation parameter. This term is called the 
Ursell parameter and has been Shown to govern the transfor- 
Mation of water waves in Shallow water. For example, it has 
been demonstrated that the nonlinear shoaling of water waves 
is predicted as the function of the Ursell's parameter. 

When Ur << 1, the linear small amplitude wave theory 
applies. In princifle, more and more terms of the fower 
series are reguired in order to keep the same relative accu- 
racy as the Ursell parameter increases. 

The Ursell parameter 1s ausefull guide, but is not 
necessarily sufficient for judging the relative importance 
of the nonlinear effects. A qualitative idea of the imper- 
tance of noniinearity is also given by the ratio of spectral 
values of each harmenic frequency to the spectral value of 
the fundamental frequency. However, the nonlinear effects of 
Shoaling transformaticn are cumulative, so that the ncnli- 
nearities may be ilccally weak, Ditwwsa Giri Cantwec ross— 
spectral transfer can occur if the Sshoaling region is wide 


in Compariscn to the nonlinear interaction distance. 
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Tift. EXPERIMENTAL PROCEDURES 

Wave and velocity data acquired at Santa fBartara, 
Califcrria, during the 2nd, 3rd, 4th, 5th and 6th Febriany 
1980, as part of the Nearshore Sediment Transport Study, are 
used to be compare with the parameterized model. The field 
data, considered in this work, include only measurements 
made py the current meters outside the breaking zone. 
Current meters CO7X, CO3X and C91X (see Figure 2) are used 
on the 2nd and 6th of February; current meters C03X, C0O1X 
and COQDX are used on the 3rd of February, and C0O1X and CODX 
are used on the 4th of February. 

For each day and each instrument the following statis- 
tics were calculated: average wave height, H, the maxinun 
wave height, Hmax, the significant wave height, H1/3, the 
root mean Syuared wave height, Hrms, and the average of the 
heights cf the 1/10 highest waves, H1/10, as well as the 
associated water depths, and the surface 


elevaticn spectrun. 


A. WAVE SPECTRAL ANALYSIS 


The surface elevation spectra were calculated frem 68 
Minute current meter records. The water particle velocity 
compcnents measured Ly the current meters were recorded ona 
Special receiver/tape recorder, digitally low-pass filtered 
and sampied at 2 samples/sec. The data were tuen high-pass 
Filtered at 9.05 hz to exclude low freguency variations. 
These current data were used to infer wave heigkts. fhe 
velocity signals were convolved using linear wave theory to 
obtain surface elevations. Tne complex Fourier spectra of 


the hcrizontal velocity components were first calculated and 
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vectorially added to give V(f). The complex surface eleva- 
tion spectrun, X(f), was calculated applying the linear wave 


theory transfer function H(f) 


X(£) = H(f£) * V(f) (354) 
where 
H(£) = (sinh kd)/{QOcosh k (d+z)} (S22) 


with z = measurement cf elevation. 

The calculated energy density spectra are to be compared 
with the parametric model. Since the parametric model 
defines all the truncated Fourier series coefficients as 
Simple analytical functions in terms of the wave height H, 
wave period T, and water depth d, it was necessary to trans- 
form the energy density spectra of the sea-surface elevation 
into a discrete (spikes) spectra aS Shown in Figure 3. 

The discrete spectra were calculated by first finding 
the fundamental freguency defined as the frequency associ- 
ated with the highest energy density peak. The harmcnics 
were defined as the integer mutiples of this fundamental 
tireguency. The bandwidths associated with the energy freaks 
were defined at the frequencies corresponding to the half 
power pcoints of the energy density values. The variance of 
the peaks, Calculated as the area of the energy density 
spectrum between half power points, describe tne height of 
the spectral spike. The analysis was carried out to the 
highest harmonic for which half power points could be 
defined. The reason for defining the Landwidths at the half 
power foints is because tne half power values correspord to 
the level at which the spectral components are independent 
of one another. It was noted that the bandwidth increases 


approximately with the number of harmonic. 
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Assubing 7) = 0, the variance of the wave surface is 


given by 


/ 2 
Vanni aa 7) at (3.3) 


The parameterized variance, or energy, can be defined from 


eguation z.15 in the discrete mode 
Z 


var (n)= K re (ey A, sinh[ak(d +7) } (3.4) 





With nh =@w2,3.. se 

The wave height (and wave amplitude) used to calculate 
the parametric energy line spectrum was the average wave 
height, H,kecause it was the wave height parameter which 
cfrered test fitting to the field data. 

Ccmparisons for all days were made between the ccnverted 
actual line spectrum and the parametric line spectrur calcu- 
lated using the inptt parameters of average wave height H, 
average water depth d, and the period corresponding tc the 
defined fundamental freguency (I = 1/£), calculated for the 


68 minute measurements. 


Be. INCIVIDUAL WAVE ANALYSIS 


A second approach was to compare individual waves tneas- 
ured in the field with the parametric model predictious. 
The £iltered and digitized velocity signals were convclved 
uSing linear theory tc obtain surface elevations. The indi- 
vidual waves were defined using the zero-up crossing method 
as expiained in Figure 10. The highest maximum and the 
lowest minimum of the surface elevation within a fperiod 
interval define the crest and the trough of a wave. A wave 
height Hi is defined as the range of f) (t) an “thag 
interval, the period is the time interval between two 


consecutive zero-up crossings cf 7) (fie 
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A Fourier analysis was performed on each individual wave 
using a Fast Fourier Transform algorithm. The parametric 
coefficients associated with the surface elevation formula 
were also determined for each wave using the measured di, 
i, and di. Statistics of Fourier coefficients from field 


data and the parametric model were generated for each day 


Nw 
En = (1/Nw) * ) An (i) (3.5) 
iz/ 
(where Nw = number of analyzed waves). For each run compar- 
isons were made between this approach and the wave spectral 
Poa lySis. Since tke average wave period was akout 15 
seconds, the total number of waves in each 68 minutes record 


was about 300, which give reasonable wave statistics. 


C. SUMMARY 


Beth analyses were based on the similarity between the 


Fourier series representation cr the sea surface elevation 


7 @) 
if (t) = a0 + > an cos ng (376) 
N=] 
(with a0 = mean) and the sea surface elevation descrited by 
the parametric model as defined by the equation 2.15. The 
coefficients to be ccmpared are the ants for the Fourier 
series defined by the Fourier integrals and tke Le Mehaute 
coefficients | 
- ({l/ JO **2) *a*g*k* An * Sinh (nk (d+7) )} 
where An are defined farametrically as function of the depth 
parameter 6 and wave steepness %6 sand 7) defined itera- 


tively from equation z.15. 
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The spectral data analysis is a linear approach. the 
spectral components are assumed independent of each other. 
The analysis considers that the spectral component fhase 
angles at each freguency are random (putting the freblenm in 
a probabilistic or stcchastic setting). On the other hand, 
the individual wave analysis is also a linear approach, kLcut 
the phase relationshipr between harmonics is presumed deter- 
Ministic (i.e. phase coupled with the fundamental frequency, 
aS evidenced by the observed peaked crests and flattened 


troughs). 
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A. SFECTRAL ANALYSIS 


1. Nonlinear Statistics 


A typical energy density spectrum for all analyzed 
instruments is presented in Figure 2 (upper panel). fhe 
presence of strong harmonics in the spectrum indicates the 
importance of nonlinearities of the waves in this region of 
Shallcw water. A qualitative idea orf the importance of 
nonlinearity is given by the ratio of the harmonic energies 
to the primary frequency energy E(nf)/E(f), where n= 2 3 
eee are the harmonic frequencies. Due to the difficulty of 
defining the half power points for the harmonics higher than 
the third, only the ratios involving tne first and seccnd 
harmonics are presented. 

Tables I - V give the ratios for ali five days 
studied for both the farameterized model (denoted L.M.) = and 
the field data (denoted by F.).The water depth and the 
Ursell parameter associated with each measurement are also 
presented. The range of water depths is from about 670 cm 
to 138 cm. The range cf the Ursell parameter is from abcut 6 
ho 1zC. 

For large Ursell numbers ( Ur > 75), the parameter- 
ized mcdel overestimates the importance of the energy asso- 
ciated with the harmonics relatively to the fundamental 
frequency energy. The same ratios for the field data are 20 
- 60 percent lower. For smali Ursell numbers ( Ur < 25), 
the parameterized model underestimates the importance of the 
harmonics relatively to the fundamental frequency energy. 
The same ratios for the field data are 0.5 - 40 fercent 


higher. 
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The best fitting between the parameterized model and 
the field data occurs at intermediate Ursell numbers ( 25 < 
Ur < 75 ). It is found that the parameterized model energy 
ratics ccmpared with the field data energy ratios have an 
error of less than 20 percent for tne cases of intermediate 
Ursell number. Some exceptions were observed to the above 
indicated percentages, but the energy contribution at these 


frequencies is very small. 


2- Local Spectral Comparisons 


An absolute ccmparison was made between the parane- 
terized model spectral energy peaks and field data spectral 
energy peaks as a function of frequency. Typical results for 
three different Ursell numbers are presented in Figures 4 
ands 25. In regions of strong nonlinear effects ( large 
Ursell number), the parameterized model underestimates the 
energy at the fundamental frequency. The energy drops 
linearly from the fundamental freyuency to the harmonics, in 
the parameterized model, while field data spectra fall cff 
much faster. In regicns of small Ursell numbers the parame- 
terized model concentrates almost all the energy in the 
fundamental frquency, overestimating it relative to the 
field data. The best agreement between the parameterized 
spectrai energy peaks and field data energy feaks are 
observed for intermediate Ursell numbers. 

The Fourier ccefficients used to evaluate field data 
spectra are plotted against the Le Mehaute spectral coeffi- 
cients ( 1/0 **2) agkA sinh {nk(d + 7))} (Figs. 6-8). A 
reasonaktle correlaticn is found for tne first three 
coefficients. 

The spectral analysis approach iS summarized by 
plotting the ratio of the Le Mehaute spectral peaks to field 
data spectral peaks against the Ursell parameter (Fig. 9). 


The ratio for the fundamental freguency energy decreases as 
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the Ursell parameter increases emphasizing how the parame- 
terized model underestimates the fundamental frequency spec- 
tral peak in situations of strong nonlinearities ( large 
Ursell parameter). MThe ratios for both the first and second 
harmonic frequency energies increase as the Ursell farameter 
increases; it must be noted, as indicated earlier, that the 


data cnly ranges in Ursell numbers from about 6 to 120. 


B. INDIVIDUAL WAVE ANALYSIS 


Individual waves were identified using the zero-up-cross 
method. A Fourier analysis was performed on each individual 
wave. The parameterized spectral coefficients were also 
determined for each individual wave. 

The Fourier coefficients associated with the fundamental 
frequency, first and second harmonics are plotted serfarately 
against the corresponding parameterized coefficients. The 
coefficients at the fundamental ‘frequency (Fig. 14) are 
reasonably well correlated (correlation coefficient = 
we o0) . The Fourier coefficients plotted against the Le 
Mehaute coefficients follow a straight line with intercept 
at the origin and slcre of abcut one. The group of plctted 
points with a greater slope (approximately 2) are associated 
with measurements in shallower waters with higher JUrsell 
numbers (see Figs. 12-16). The graphs of the the first and 
second harmonic coefficients ({Figs. 17-18) show similar 
Characteristics as the fundamental frequency; however, the 
correlation between the coefficients is poorer. 

The ratio of the Ile Mehaute coefficients to the Fourier 
coefficients is alsc plotted against the JUrsell parameter 
The ratio for the fundamental frequency coefficients 
converge to a bound vaiue of 0.5 at larye Ursell numbers 
(Fig. 19). For Ursell numbers less than 75 the ratio of the 


Le Mehaute to Fourler coefiicients is akout 1 with the 
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exception of a few pcints near Ur = 0 which were identified 
to be related with noisy Signals. For large Ursell nunrbers 
the Fourier coefficients are about twice the Le Mehaute 
coefficients, again showing the tendency that the parameter- 
ized model has to underestimate the fundamental peak energy 
in regions of strong nonlinearities (large Ursell numbers). 
The ratios for the first and second harmonic coeffi- 
cients (Figs. 19 and 20) show a bounded value near unity as 
the ursell parameter increases. Noisy signals near Ur = 0 


are also cbserved. 


C. CCMPARISONS BETWEEN APPROACHES 


Both linear (spectral analysis) and nonlinear (indi- 
vidual wave analysis) approaches lead to Similar results. 
The parameterized mcdel is observed to overestimate the 
fundamental frequency energy in regions of small JUrsell 
number (weak nonlinearities) and underestimate it in regions 
of large Ursell number. The individual wave analysis fres- 
ents a limit tendency for large JUrsell numbers, with the 
fundamental frequency Fourier coefficients about twice those 
of the parameterized spectral coefficients. The better 
correlaticn observed between the Fourier coefficients and 
the Le Mehaute coefficients in the individual wave analysis 
as compared with spectral analysis approach is due, ferhars, 
to the definition of half-power criteria or to the averaging 
processes involved in calculating the wave spectra. 

The reason for the siMilarity of results for the two 
approaches can be seen by comparing the spectra caiculated 
in the usual manner with a "ncenlinear spectrun". A nenli- 
near spectrum was calculated uSing the Fourier amplitude 
coefficients calculated for the individual waves. Using the 
same Fandwidths as the ordinary spectrum, amplitude coeffi- 
cients were sorted into the respective freguency bands anda 


nonlinear srectrum calculated as 


az 


Ny 
Gnl (£) = (1/Nw) ) a? (£) VON (4.1) 
i=/ 
where the coefficients Ai(f) are centered on frequency f 
Within bandwidth Af. The total number of waves is denoted 
by Nw. An example of the ccmfarison of the nonlinear spec- 
trum calculated in this manner with the linear spectral 
method is shown in Figure 21. The two spectra are Similar 
and bcth show strong harmonic peaks. The spectral energy in 
the ordinary spectrum tends to be smeared more, whereas the 
peaks and valleys in the nonlinear spectrum are more sharfly 
defined. 
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¥V. CONCLUSIONS 





The coefficients of the fparameterized solution tc the 
monochromatic, irrotational, nonlinear wave over a hori- 
zontal botton, as developed by Le Mehaute et al. (1984), are 
expressed in terms cf simple analytical functions of wave 
height, H, wave period, T,and water depth, d. All gquanti- 
ties ~-e€xpiaei tly derived fron the potential function 
(velocity components, dynamic fressure, acceleration, i1.¢é., 
the internal wave field) can be directly determined. The 
relative simplicity cf the calculations of the parameterized 


soluticn recommends it for engineering applications. 


The wave spectra determined from the parameterized solu- 
tion show good agreement with the measured wave spectra from 
field data for conditions of intermediate Ursell numkers (25 
< Ur < 75). For small Ursell numbers, the parameterized 
soluticn concentrates almost all the energy in the funda- 
mental frequency representing a sinusoidal waveforn. Under 
these conditions, the field wave spectra always contained 
relatively less energy associatead with the fundarental 


frequency and more energy distributed to the harmonics. 


Fer large Ursell numbers, the parameterized solution 
descrites a waveform with more peaked crests and flatter 
throughs relative to the measured waves, Since it underesti- 
mates the energy asscciated with the fundamental frequency 
and overestimates tke energy associated with the higher 


narmonics. 


The ratios between the Le Mehaute and fieid data speéec- 
tral coefficients associated with the ist and 2nd harmonics 


present a bound ligit tendency of 1 for large Urseil 
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humbers, indicating gcod agreement between the two processes 
for these conditions. The same ratios for the coefficients 
associated with the the fundamental frequency show a Lound 
limit tendency of 0.5, indicating that the first coeffi- 
cient of the parameterized model needs to be corrected for 


strong (Ur > 75) nonlinear effects. 


Further compariscns should be made for different sea 
conditions to determine the corrections to the parameterized 
model coefficients. However, it must be noted that a liri- 
tation of the parameterized solution is that its accuracy 
decreases rapidly as the wave heights approach limiting wave 


conditions due to the use of a truncated Fourier series. 
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Figure 1 Flow Chart of the Parameterized Nonlinear 
Wave Model. 
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SANTA BARBARA BEACH 
CALIFORNIA 
2-6 FEB. 1980 


WAVE STAFFS 





Figure 2 Beach Profile with the Instruments Location. 
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Figure 3 Typacat fue Density Spectrum (above) and 
Ccnversion into [Liscrete Energy Spectrum (below). 
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TABLE I 


Ratio of Harmonic Ener to Primary Frequency Energy 
037 PEE 80 , 


DEPTH (cm) idle 308 Sholay bie log levees 
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TABLE It 


Ratio of Harmonic Ener to Primary Frequency Ener 
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TABLE I1ILf 


atio of Harmonic Ener to Primary Frequency Ener 
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DEPTH (en 270.8 
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TABLE IV 
Ratio of Harmonic Ener to Primary Frequency Ener 
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TABLE V 


Ratio of Harmonic Ener to Primary Frequency Ener 
OG FER. 80 Y 3 ¥ 7 


DEPTH (cm) 138.8 2B 7 345.4 
eee Ee 


iam Dor | a a ce | 
cr )/ECE) 0.548 Ole, eis ial U2 196 O02 we 110 
icon )/ E(t) Deze 0.068 02009 0.088 He002 0.049 
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Figure 5 local Spectral Comparisons for 
Intermediate Ursell Number. 
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where X 2S a Crest 


~ is a trough 


O is the zero-up crossing 


ie lee ina a inn i a ch le a Ga 
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Figure 10 Defining Waves using 
Zero-up Cross Method. 
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